Comparison of parametric, orthogonal, and spline functions to model individual lactation curves for milk yield in Canadian Holsteins

Test day records for milk yield of 57,390 first lactation Canadian Holsteins were analyzed with a linear model that included the fixed effects of herd-test date and days in milk (DIM) interval nested within age and calving season. Residuals from this model were analyzed as a new variable and fitted with a five parameter model, fourth-order Legendre polynomials, with linear, quadratic and cubic spline models with three knots. The fit of the models was rather poor, with about 30%-40% of the curves showing an adjusted R-square lower than 0.20 across all models. Results underline a great difficulty in modelling individual deviations around the mean curve for milk yield. However, the Ali and Schaeffer (5 parameter) model and the fourth-order Legendre polynomials were able to detect two basic shapes of individual deviations among the mean curve. Quadratic and, especially, cubic spline functions had better fitting performances but a poor predictive ability due to their great flexibility that results in an abrupt change of the estimated curve when data are missing. Parametric and orthogonal polynomials seem to be robust and affordable under this standpoint

Non-additive genetic effects for fertility traits in Canadian Holstein cattle (Open Access publication )

The effects of additive, dominance, additive by dominance, additive by additive and dominance by dominance genetic effects on age at first service, non-return rates and interval from calving to first service were estimated. Practical considerations of computing additive and dominance relationships using the genomic relationship matrix are discussed. The final strategy utilized several groups of 1000 animals (heifers or cows) in which all animals had a non-zero dominance relationship with at least one other animal in the group. Direct inversion of relationship matrices was possible within the 1000 animal subsets. Estimates of variances were obtained using Bayesian methodology via Gibbs sampling. Estimated non-additive genetic variances were generally as large as or larger than the additive genetic variance in most cases, except for non-return rates and interval from calving to first service for cows. Non-additive genetic effects appear to be of sizeable magnitude for fertility traits and should be included in models intended for estimating additive genetic merit. However, computing additive and dominance relationships for all possible pairs of individuals is very time consuming in populations of more than 200 000 animals

Fit of different functions to the individual deviations in random regression test day models for milk yield in dairy cattle

The shape of individual deviations of milk yield for dairy cattle from the fixed part of a random regression test day model (RRTDM) was investigated. Data were 53,217 TD records for milk yield of 6,229 first lactation Canadian Holsteins in Ontario. Data were fitted with a model that included the fixed effects of herd-testdate, DIM interval nested within age and season of calving. Residuals of the model were then fitted with the following functions: Ali and Schaeffer 5 parameter model, fourth-order Legendre Polynomials, and cubic spline with three, four or five knots. Result confirm the great variability of shape that can be found when individual lactation are modeled. Cubic splines gave better fitting pe4rformances although together with a marked tendency to yield aberrant estimates at the edge of the lactation trajectory

Identification of unique ROH regions with unfavorable effects on production and fertility traits in Canadian Holsteins

Background: The advent of genomic information and the reduction in the cost of genotyping have led to the use of genomic information to estimate genomic inbreeding as an alternative to pedigree inbreeding. Using genomic measures, effects of genomic inbreeding on production and fertility traits have been observed. However, there have been limited studies on the specific genomic regions causing the observed negative association with the trait of interest. Our aim was to identify unique run of homozygosity (ROH) genotypes present within a given genomic window that display negative associations with production and fertility traits and to quantify the effects of these identified ROH genotypes. Methods: In total, 50,575 genotypes based on a 50K single nucleotide polymorphism (SNP) array and 259,871 pedigree records were available. Of these 50,575 genotypes, 46,430 cows with phenotypic records for production and fertility traits and having a first calving date between 2008 and 2018 were available. Unique ROH genotypes identified using a sliding-window approach were fitted into an animal mixed model as fixed effects to determine their effect on production and fertility traits. Results: In total, 133 and 34 unique ROH genotypes with unfavorable effects were identified for production and fertility traits, respectively, at a 1% genome-wise false discovery rate. Most of these ROH regions were located on bovine chromosomes 8, 13, 14 and 19 for both production and fertility traits. For production traits, the average of all the unfavorably identified unique ROH genotypes effects were estimated to decrease milk yield by 247.30 kg, fat yield by 11.46 kg and protein yield by 8.11 kg. Similarly, for fertility traits, an average 4.81-day extension in first service to conception, a 0.16 increase in number of services, and a - 0.07 incidence in 56-day non-return rate were observed. Furthermore, a ROH region located on bovine chromosome 19 was identified that, when homozygous, had a negative effect on production traits. Signatures of selection proximate to this region have implicated GH1 as a potential candidate gene, which encodes the growth hormone that binds the growth hormone receptor. This observed negative effect could be a consequence of unfavorable alleles in linkage disequilibrium with favorable alleles. Conclusions: ROH genotypes with unfavorable effects on production and fertility traits were identified within and across multiple traits on most chromosomes. These identified ROH genotypes could be included in mate selection programs to minimize their frequency in future generations

Parent-offspring genotyped trios unravelling genomic regions with gametic and genotypic epistatic transmission bias on the cattle genome

Several biological mechanisms affecting the sperm and ova fertility and viability at developmental stages of the reproductive cycle resulted in observable transmission ratio distortion (i.e., deviation from Mendelian expectations). Gene-by-gene interactions (or epistasis) could also potentially cause specific transmission ratio distortion patterns at different loci as unfavorable allelic combinations are under-represented, exhibiting deviation from Mendelian proportions. Here, we aimed to detect pairs of loci with epistatic transmission ratio distortion using 283,817 parent-offspring genotyped trios (sire-dam-offspring) of Holstein cattle. Allelic and genotypic parameterization for epistatic transmission ratio distortion were developed and implemented to scan the whole genome. Different epistatic transmission ratio distortion patterns were observed. Using genotypic models, 7, 19 and 6 pairs of genomic regions were found with decisive evidence with additive-by-additive, additive-by-dominance/dominance-by-additive and dominance-by-dominance effects, respectively. Using the allelic transmission ratio distortion model, more insight was gained in understanding the penetrance of single-locus distortions, revealing 17 pairs of SNPs. Scanning for the depletion of individuals carrying pairs of homozygous genotypes for unlinked loci, revealed 56 pairs of SNPs with recessive epistatic transmission ratio distortion patterns. The maximum number of expected homozygous offspring, with none of them observed, was 23. Finally, in this study, we identified candidate genomic regions harboring epistatic interactions with potential biological implications in economically important traits, such as reproduction

The genetic architecture of milk ELISA scores as an indicator of Johne's disease (paratuberculosis) in dairy cattle

peer-reviewedJohne's disease (or paratuberculosis), caused by Mycobacterium avium ssp. paratuberculosis (MAP) infection, is a globally prevalent disease with severe economic and welfare implications. With no effective treatment available, understanding the role of genetics influencing host infection status is essential to develop selection strategies to breed for increased resistance to MAP infection. The main objectives of this study were to estimate genetic parameters for the MAP-specific antibody response using milk ELISA scores in Canadian Holstein cattle as an indicator of resistance to Johne's disease, and to unravel genomic regions and candidate genes significantly associated with MAP infection. After data editing, 168,987 milk ELISA records from 2,306 herds, obtained from CanWest Dairy Herd Improvement, were used for further analyses. Variance and heritability estimates for MAP infection status were determined using univariate linear animal models under 3 scenarios: (a) SCEN1: the complete data set (all herds); (b) SCEN2: herds with at least one suspect or test-positive animal (ELISA optical density ≥0.07); and (c) SCEN3: herds with at least one test-positive animal (ELISA optical density ≥0.11). Heritability estimates were calculated as 0.066, 0.064, and 0.063 for SCEN1, SCEN2, and SCEN3, respectively. The correlations between estimated breeding values for resistance to MAP infection and other economically important traits, when significant, were favorable and of low magnitude. Genome-wide association analyses identified important genomic regions on Bos taurus autosome (BTA)1, BTA7, BTA9, BTA14, BTA15, BTA17, BTA19, and BTA25 showing significant association with MAP infection status. These regions included 2 single nucleotide polymorphisms located 2 kb upstream of positional candidate genes CD86 and WNT9B, which play key roles in host immune response and tissue homeostasis. This study revealed the genetic architecture of MAP infection in Canadian Holstein cattle as measured by milk ELISA scores by estimating genetic parameters along with the identification of genomic regions potentially influencing MAP infection status. These findings will be of significant value toward implementing genetic and genomic evaluations for resistance to MAP infection in Holstein cattle

Association of genetic polymorphisms related to Johne’s disease with estimated breeding values of Holstein sires for milk ELISA test scores

peer-reviewedBackground Johne’s disease (JD) is a chronic intestinal inflammatory disease caused by Mycobacterium avium subsp. paratuberculosis (MAP) infection in ruminants. Since there are currently no effective vaccine or treatment options available to control JD, genetic selection may be an alternative strategy to enhance JD resistance. Numerous Single Nucleotide Polymorphisms (SNPs) have been reported to be associated with MAP infection status based on published genome-wide association and candidate gene studies. The main objective of this study was to validate these SNPs that were previously identified to be associated with JD by testing their effect on Holstein bulls’ estimated breeding values (EBVs) for milk ELISA test scores, an indirect indicator of MAP infection status in cattle. Results Three SNPs, rs41810662, rs41617133 and rs110225854, located on Bos taurus autosomes (BTA) 16, 23 and 26, respectively, were confirmed as significantly associated with Holstein bulls’ EBVs for milk ELISA test score (FDR < 0.01) based on General Quasi Likelihood Scoring analysis (GQLS) analysis. Single-SNP regression analysis identified four SNPs that were associated with sire EBVs (FDR < 0.05). This includes two SNPs that were common with GQLS (rs41810662 and rs41617133), with the other two SNPs being rs110494981 and rs136182707, located on BTA9 and BTA16, respectively. Conclusions The findings of this study validate the association of SNPs with JD MAP infection status and highlight the need to further investigate the genomic regions harboring these SNPs

Reproductive tract size and position score: Estimation of genetic parameters for a novel fertility trait in dairy cows.

The dairy industry is moving toward selecting animals with better fertility to decrease the economic losses linked to reproductive issues. The reproductive tract size and position score (SPS) was recently developed in physiological studies as an indicator of pregnancy rate and the number of services to conception. Cows are scored as SPS 1, 2, or 3 based on the size of their reproductive tract and its position in the pelvis, as determined by transrectal palpation. The objective of this study was to estimate genetic parameters for SPS to assess its potential as a novel fertility trait. Phenotypes were collected at the University of British Columbia's research herd from 2017 to 2020, consisting of 3,247 within- and across-lactation SPS records from 490 Holstein cows. A univariate animal model was used to estimate the variance components for SPS. Both threshold and linear models were fit under a Bayesian approach and the results were compared using the Spearman rank correlation (r) between the estimated breeding values. The 2 models ranked the animals very similarly (r = 0.99), and the linear model was selected for further analysis. Genetic correlations with other currently evaluated traits were estimated using a bivariate animal model. The posterior means (± posterior standard deviation) for heritability and repeatability within- and across-lactation were 0.113 (± 0.013), 0.242 (± 0.012), and 0.134 (± 0.014), respectively. The SPS showed null correlations with production traits and favorable correlations with traditional fertility traits, varying from -0.730 (nonreturn rate) to 0.931 (number of services). Although preliminary, these results are encouraging because SPS seems to be more heritable than and strongly genetically correlated with number of services, nonreturn rate, and first service to conception, indicating potential for effective indirect selection response on these traits from SPS genetic selection. Therefore, further studies with larger data sets to validate these findings are warranted

Predicting dry matter intake in Canadian Holstein dairy cattle using milk mid-infrared reflectance spectroscopy and other commonly available predictors via artificial neural networks.

Dry matter intake (DMI) is a fundamental component of the animal's feed efficiency, but measuring DMI of individual cows is expensive. Mid-infrared reflectance spectroscopy (MIRS) on milk samples could be an inexpensive alternative to predict DMI. The objectives of this study were (1) to assess if milk MIRS data could improve DMI predictions of Canadian Holstein cows using artificial neural networks (ANN); (2) to investigate the ability of different ANN architectures to predict unobserved DMI; and (3) to validate the robustness of developed prediction models. A total of 7,398 milk samples from 509 dairy cows distributed over Canada, Denmark, and the United States were analyzed. Data from Denmark and the United States were used to increase the training data size and variability to improve the generalization of the prediction models over the lactation. For each milk spectra record, the corresponding weekly average DMI (kg/d), test-day milk yield (MY, kg/d), fat yield (FY, g/d), and protein yield (PY, g/d), metabolic body weight (MBW), age at calving, year of calving, season of calving, days in milk, lactation number, country, and herd were available. The weekly average DMI was predicted with various ANN architectures using 7 predictor sets, which were created by different combinations MY, FY, PY, MBW, and MIRS data. All predictor sets also included age of calving and days in milk. In addition, the classification effects of season of calving, country, and lactation number were included in all models. The explored ANN architectures consisted of 3 training algorithms (Bayesian regularization, Levenberg-Marquardt, and scaled conjugate gradient), 2 types of activation functions (hyperbolic tangent and linear), and from 1 to 10 neurons in hidden layers). In addition, partial least squares regression was also applied to predict the DMI. Models were compared using cross-validation based on leaving out 10% of records (validation A) and leaving out 10% of cows (validation B). Superior fitting statistics of models comprising MIRS information compared with the models fitting milk, fat and protein yields suggest that other unknown milk components may help explain variation in weekly average DMI. For instance, using MY, FY, PY, and MBW as predictor variables produced a predictive accuracy (r) ranging from 0.510 to 0.652 across ANN models and validation sets. Using MIRS together with MY, FY, PY, and MBW as predictors resulted in improved fitting (r = 0.679-0.777). Including MIRS data improved the weekly average DMI prediction of Canadian Holstein cows, but it seems that MIRS predicts DMI mostly through its association with milk production traits and its utility to estimate a measure of feed efficiency that accounts for the level of production, such as residual feed intake, might be limited and needs further investigation. The better predictive ability of nonlinear ANN compared with linear ANN and partial least squares regression indicated possible nonlinear relationships between weekly average DMI and the predictor variables. In general, ANN using Bayesian regularization and scaled conjugate gradient training algorithms yielded slightly better weekly average DMI predictions compared with ANN using the Levenberg-Marquardt training algorithm